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Abstract 

We present some recent advances in the computational study of baryon num- 
ber violation in high energy electroweak collisions. We examine classically allowed 
processes above the sphaleron barrier and, using a stochastic search procedure, we 
explore the topology changing region in the energy and particle number plane. Find- 
ing topology changing classical solutions with small incident particle number would 
be an indication that baryon number violation becomes unsuppressed in high energy 
collisions. Starting with a topology changing solution of approximately 50 incoming 
particles, our Monte-Carlo procedure has produced other topology changing solutions 
with 40% lower incident particle numbers, with energies up to one and a half times the 
sphaleron energy. While these solutions still involve a rather large number of incident 
particles, we have nonetheless demonstrated that our search procedure is effective in 
reducing the particle number while ensuring topology change. Taking advantage of 
more powerful computational resources, we plan to extend the search to still higher 
energies. 
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1 Introduction 



The prospect of observable high energy baryon number violation within the stan- 
dard model has recently attracted widespread attention. Unfortunately, despite con- 
siderable effort by a great many theorists, the issue still remains largely unsettled. 
The purpose of this paper is to explain some recent computational developments that 
shed more light on the problem and which might help contribute to a final solution. 
Before we begin, however, in an effort to write a self-contained work, we shall give a 
brief exposition of nonperturbative baryon number violation in the standard model. 
We unfortunately cannot survey the vast literature on the subject with the depth it 
deserves, so instead a brief summary of the facts germane to our numerical approach 
must suffice. 

For our purposes, when we talk of the "standard model" we mean the standard 
model in which the Weinberg angle has been set to zero, i.e. we shall be considering 
SU{2) gauge theory spontaneously broken via a single Higgs doublet. This simplified 
model has all the relevant physics. Most importantly, the gauge structure dictates 
nontrivial topology for the bosonic vacuum sector. Working in the temporal gauge 
with periodic boundary conditions at spatial infinity, each vacuum may be character- 
ized by an integer called the winding number which measures the number of times 
the gauge manifold is wound around 3-space]l|. As this number is a topological in- 
variant, vacua of different winding numbers cannot be continuously deformed into 
one another. 

Because of the axial vector anomaly, baryon number violation occurs when the 
gauge and Higgs fields change their topology p[. Different topological sectors are 
separated by an extremely high barrier, the top of which is a static saddle-point 
solution to the equations of motion. This configuration is called the sphaleron0, and 
it has an energy of about 10 TeV and possesses a single unstable direction in field 
space. At low energy the baryon number violating rates are exceedingly small, as the 
gauge and Higgs fields must first tunnel through the sphaleron, which is extremely 
unlikely indeed. 

Recently the prospect of rapid baryon number violation at high temperatures 
and high energies has emerged. The basic idea is that if the gauge and Higgs fields 
have enough energy, they can change their topology by sailing over the sphaleron 
barrier rather than being forced to tunnel through it. At high temperatures this 
is precisely what happens, and it is generally agreed that baryon number violation 
becomes unsuppressed in the early universe @]. 

The situation in high energy collisions is far less clear. The limiting process in 
baryon number violation is the production of a sphaleron-like configuration from an 
incident beam of high energy particles. But since the sphaleron is a large extended 
object, there is a scale mismatch with the initial high energy 2-particle state, and 
hence one naively expects the baryon number violating rate to be small. However, 



Ringwald and EspinosaQ have suggested that the sum over many-particle final 
states gives rise to factors that grow with energy sufficiently rapidly to offset any 
exponential suppression. If true, this offers the exciting prospect of one day studying 
baryon number violation in the laboratory. Their approach, however, neglects some 
important corrections which still elude calculation despite considerable effort. Apart 
from lattice simulations, semi-classical techniques are our only handle on nonpertur- 
bative effects. The basic problem with anomalous baryon number violation is that 
exclusive two-particle initial states are not very amenable to these methods, and there 
are potentially large initial state corrections whose effects remained undetermined. 

Our efforts lie in an attempt to overcome these limitations. Many people have 
struggled in similar endeavors, but we only have space to summarize the work of one 
group upon which our approach has been partly inspired. In an effort to alleviate 
problems arising from exclusive two-particle states, Rubakov, Son and TinyakovQ 
consider an inclusive quantity: 

a{E,N) = \<f\SPEPNV>? (1) 

where the sum is over all initial and final states, S is the S'-matrix, and Pe and P/v are 
projection operators onto subspaces of energy E and particle number N respectively. 

Unlike the exclusive two-particle cross section (J2{E), the quantity a{E,N) is 
directly calculable by semiclassical methods. If the energy and particle number are 
parameterized by 

E = (2) 
9 

N = (3) 
9^ 

then in the limit g —>■ with e and u held fixed, the path integral for a{E, N) can be 
saturated by a classical saddle-point solution to the equations of motion. The cross 
section takes the form 



a{E, N) = exp 



lF{e,u) + Oig') 
9^ 



(4) 



where the function F{e, v) is determined by the classical solution. These solutions 
naturally divide into two regimes: there is Euclidean evolution corresponding to tun- 
neling under the sphaleron barrier, and Minkowski evolution corresponding to classical 
motion before and after the tunneling event. 

The utility of a{E, N) arises because it may be used to bound a2{E). By construc- 
tion, a{E,N) provides an upper bound to a2{E). A lower bound may be obtained 
under some reasonable physical assumptions. The two-particle process is expected to 
be dominated by a most favored intermediate sphaleron-like state, and the rate into 



this intermediate state bounds the two-particle cross section from below. Combining 
these upper and lower bounds allows one to write |0, || 



exp{-const N) (t{E, N) < a2{E) < (r{E, N) , (5) 
from which it follows that[0] 

lima2(^) = F(e,z/) + 0(z/) . (6) 

The consistency of the first inequality requires that F{e, u) have a smooth z/ ^ 
limit, in which case F{e,0) determines the exponential behavior of a2{E). However, 
the second inequality of (H) holds regardless of continuity, and hence if a{E,N) is 
exponentially suppressed then so is o"2(-E). 

Obviously finding these Euclidean-Minkowski hybrid solutions will be extremely 
illuminating, and we are presently engaged in this rather formidable numerical task. 
In this paper, however, rather than exploring the full barrier penetration problem, 
we report on a complimentary approach. We examine the classically allowed regime 
above the sphaleron barrier in which the saddle-points that saturate the path integral 
are pure Minkowski solutions. This is less computationally demanding than solving 
the tunneling problem, while still yielding considerable information about baryon 
number violation. Spatial limitations prevent us from giving a full blown treatment 
of our numerical investigation, and the reader is referred to Ref. |jlO|] for complete 



details. But the basic idea is that if a topology changing classical solution with 
small incident particle number could be found, this would be a strong indication that 
baryon number violation would be observable in high energy two-particle collisions. 
Conversely, if there are no small-multiplicity topology changing solutions, then it is 
unlikely that the rates become exponentially unsuppressed. 

This can be made more precise in the following manner. Because of energy dissi- 
pation, the system will asymptotically approach vacuum values and will consequently 
linearize in the past and future. Field evolution then becomes a superposition of 
normal mode oscillators with amplitudes a„, which allows us to define the asymptotic 
particle number i/ = J^Wnl"^- Furthermore, since the fields approach vacuum values 
in the infinite past and future, the winding numbers of the asymptotic field config- 
urations are also well defined, and fermion number violation is given by the change 
in topology of these vacua[|ll|. Because of the sphaleron barrier, classical solutions 



that change topology must have energy e greater than that of the sphaleron. The 
problem we would like to solve, then, is whether the incident particle number u of 
these solutions can be made arbitrarily small. That is to say, we wish to map the 
region of topology changing classical solutions in the e-u plane. 

We could easily parameterize incoming configurations in terms of small pertur- 
bations about a given vacuum, but it would be extremely difficult to choose the 
parameters to ensure a subsequent change in winding number. This is because topol- 
ogy changing classical solutions must pass over the sphaleron barrier at some point 



in their evolution, which is extremely difficult arrange by an appropriate choice of 
initial conditions. So computationally we purse a different strategy. We shall evolve 
a configuration near the top of the sphaleron barrier until it linearizes and the particle 
number can be extracted. The time reversed solution, then, has a known incident 
particle number and will typically pass over the sphaleron barrier thereby changing 
topology. Of course we have no obvious control over the asymptotic particle number 
of the initial sphaleron- like configuration; however, by using suitable stochastic sam- 
pling techniques, we can systematically lower the particle number while ensuring a 
change of topology. This will allow us to explore the e-v plane and map the region 
of topology change, the lower boundary of which should tell us a great deal about 
baryon number violation in high energy collisions. 



2 Topological Transitions 



Let us now put some flesh on the bones of the introductory discussion. For 
simplicity we consider the standard model with the Weinberg angle set to zero. The 
resulting spontaneously broken SU (2) gauge theory possesses all the relevant physics 
while undergoing notable simplification. The action for the bosonic sector of this 
theory is 

S = Jdx^ j-iTrF^^F'^^ + D^^tL"^^ - A($t$ - 1)^1 , (7) 

where the indices run from to 3 and where 

= d^A,-d,A^-i[A^,A,] (8) 
D^^ = {d^-iA^)^. (9) 

We use the standard metric t] = diag(l, — 1, — 1, — 1), and have eliminated several 
constants by a suitable choice of units. We have also set (? = 1, but we shall restore 
the gauge coupling to its physical value of g = 0.652 when needed. For our numerical 
work we take A = 0.1, which corresponds to a Higgs mass of about = 72 GeV. 

To yield a computationally manageable system, we work in the spherical Ansatz 
of Ref. |12[ in which the gauge and Higgs fields are parameterized in terms of six real 



functions , ai , a , P , fi and u oi r and t: 

Ao(x,t) = ^ao{r,t)(T-it (10) 

Ai(x,t) = - [ai(r, t)(T ■ + ^^(a' - cr ■ M') + ^ ^ ^K '^'^x^a''] (11) 
2 r r 

<l>(x,t) = [/i(r,t)+Mr,t)cr-x]e , (12) 



where x is the unit three-vector in the radial direction and ^ is an arbitrary two- 
component complex unit vector. For the 4-dimensional fields to be regular at the 



origin, ao, a, ai — a/r, (1 + /3)/r and v must vanish like some appropriate power of 
r as r ^ 0. 

These spherical configurations reduce the system to an effective 1+1 dimensional 
theory whose action can be obtained by inserting (p!0|)-([I^) into from which one 
obtains f|T2| 



S = 4:7r I dt J dr 



{\X?- if - \{\X? + m? - Re(zx*0^) (13) 



2r2 

\ 2 / I i |2 



where the indices now run from to 1 and are raised and lowered with ?7^jy 
diag(l, —1), and where 



fflU 

X 



D,.6 



a + ip 
fi + iu 
{df, -ia^)x 

{df, - ^ a^)(f) 



(14) 
(15) 
(16) 
(17) 

(18) 



(^fi vac 


= d^Q 






Xvac 


= -ze^^ 


0vac 





This is an effective 2-dimensional U{1) gauge theory spontaneously broken by two 
scalar fields. It possesses the same rich topological structure as the full 4-dimensional 
theory and provides an excellent testing ground for numerical exploration. 

Vacuum states of the effective 2-dimensional theory are characterized by |x| = 
|0| = 1 and ix* 4> = —I (as well as D^x = -0/^0 = 0). The vacua then take the form 

(19) 
(20) 
(21) 

where the gauge function Q = Q{r, t) is required to vanish at r = to ensure regularity 
of the 4-dimensional fields. Furthermore, like the full 4-dimensional theory, these 
vacua still possesses nontrivial topological structure. Compactification of 3-space 
requires that Q{r,t) 27m as r — > oo, in which case the winding number of such 
vacua in the the oq = gauge is simply the integer n. Note that as r varies from zero 
to infinity, x winds n times around the unit circle while only winds by half that 
amount. 

Since the winding number is a topological invariant, a continuous path connecting 
two inequivalent vacua must at some point leave the manifold of vacuum configura- 
tions. Along this path there will be a configuration of maximal energy, and of all such 
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Figure 1: The X fields for a vacuum-to- vacuum topology changing transition in a gauge 
inconsistent with compactified 3-space. The scalar fields are traced in the complex plane as the the 
spatial coordinate spans the entire axis. Figs, (a) and (c) represent two inequivalent topological 
vacua while (b) is the sphaleron barrier separating them. 



maximal energy configurations there exists a unique one of minimal energy 0. Tliis 
configuration is called the sphaleron and may conveniently be parameterized by 

(^sphir) = , 
Xsph(r) = ^[2/(r)-l] 
4>sph{r) = ih{r) 

where the profile functions / and h vanish at r = 0, tend to unity as r ^ oo and 
are otherwise determined by minimizing the energy functional. The energy of the 
sphaleron depends very weakly on the Higgs mass and is about M^j ~ 10 TeV, or 
e = 47r (2.54) in the units we are using. 

While this form of the sphaleron in which vanishes and x and are pure 
imaginary is convenient for numerical work, it does have a slight peculiarity that we 
briefly mention. Recall that compactification of 3-space requires the gauge function V 
to approach an even multiple of 27r as r — oo. It is possible to relax this restriction, 
and it will often be convenient to choose a gauge in which f/ — (2n-|-l)7r as r — cx3, in 
which case Xvac — > i and 0vac This is precisely the large-r boundary condition of 

the sphaleron, which illustrates that (^) is inconsistent with spatial compactification. 
There is of course nothing wrong with this, and a topological transition of unit winding 
number change in this gauge is illustrated in Fig. 1. Rather than x winding once 
around the unit circle, it instead winds over the left hemisphere before the transition 
and over the right after the transition. The total phase change is still 27r, as it must 
be since this is a gauge invariant quantity. 



(22) 



Throughout most of this paper we shall use a gauge consistent with ( p2D in which 
space cannot be compactified. From a computational perspective, this will make 
perturbations about the sphaleron more easily parameterized. There will, however, 
be times in which it is more convenient to impose spatial compactification, but we 
will always alert the reader to such a change of gauge. 



3 Classical Evolution 



So far we have primarily been considering topology changing sequences of con- 
figurations, not necessarily solutions of the equations of motion. Now we turn to 
the classical evolution of the system. We will consider solutions that linearize in 
the distant past and future, and hence ones that asymptote to specific topological 
vacua. This allows us to define the incident particle number, and it makes clear what 
is meant by topology change of a classical solution (namely, the change in winding 
number of the asymptotic vacua). 

The field equations are coupled nonlinear particle differential equations, and must 
be solved computationally on the lattice. But before we present our numerical pro- 
cedure, we first formulate the problem in the continuum. The equations of motion 
resulting from the action (|T3) are 



X 



+ 1] 



2\r' 



1 



(23) 

(24) 
(25) 



The = equation in (|2^) is not dynamical but is simply the Gauss's law constraint. 

To solve these equations, we must supplement them with boundary conditions. 
The conditions at r = can be derived by requiring the 4-dimensional fields to be 
regular at the origin. The behavior as r ^ must be 



ao 


= ao,ir + ao,3r^ + . . . 


(26) 




= aifi + ai,2?^^ + • • • 


(27) 


a 


= air + a^r^ + . . . 


(28) 


/5 


= -l + /52r2 + ... 


(29) 


/i 


= /io + fi2r^ + ■■■ 


(30) 




= vir + z/ar^ + . . . , 


(31) 



where the coefficients of the r-expansion are undetermined functions of time. The 
r-behavior of the various terms are determined by the requirement that it has the 



appropriate power of r = (x^ + + 2;^)^/^ to render the 4- dimensional fields analytic 
in terms of x, y and z. For example, Oq must be odd in r since Aq is proportional to 
aocr ■ X = {ao/r)cr ■ x. In terms of x and the boundary conditions at r = become 



ao(0,t) = (32) 

X(0,t) = (33) 

Re (9^0(0, t) = (34) 

Im0(O,t) = 0. (35) 

There is another r = boundary condition which arises from the requirement that 
ai — a/r be regular as r — 0. This condition can be written as ai^ = ai, and once 
imposed on initial configurations it remains satisfied at subsequent times because of 
Gauss's law. 

We turn now to the large-r boundary conditions. Finite energy configurations 
must approach pure vacuum at spatial infinity, and we may choose a gauge in which 

a^{r,t) ^ (36) 
X(r,t) ^ t (37) 
0(r,t) ^ I (38) 

as r — > 00. This choice of gauge does not admit spatial compactification, but nonethe- 
less it is numerically conveniently since it is consistent with the simple parameteri- 
zation of the sphaleron (P^. At times we will choose a gauge consistent with spatial 
compactification in which x(r, t) — > —i and (f){r,t) — 1 as r — 00, but unless other- 
wise specified we will take the large-r boundary conditions to be (^)-(^). 

The field equations (p3D-(p^), together with boundary conditions (|3^ ) -(l38| ) , may 



now be used to evolve initial profiles and investigate their subsequent topology change. 
The evolution is performed by discretizing the system using the methods of lattice 
gauge theory, in which we subdivide the r-axis into N equal intervals of length Ar 
with finite extent L = NAr (in our numerical simulations we take N = 2239 and 
Ar = 0.04). The field theoretic system then becomes finite and may be solved using 
standard numerical techniques. 

The fields xi'f'yt) 0(^?'^) become discrete variables Xi{t) (ptit) associated 
with the lattice sites r^ = iAr where i = ■ ■ ■ N. The continuum boundary condi- 
tions render the variables at the spatial end-points nondynamical, taking the values 
Xo = —i, Xn = and (j)^ = i (the value of 0o will be discussed momentarily). The 
time component of the gauge field ao(r, t) is also associated with the lattice sites and 
is represented by the variables a^iit) with i = - ■ ■ N. We will usually work in the 
temporal gauge in which ao.i = 0, and we will not concern ourselves with this degree 
of freedom. 

The spatial components of the gauge field ai (r, t) become discrete variables asso- 
ciated with the oriented links of the lattice, and we represent them by ai,i(t) = ai(t) 



located at positions rj_|_i/2 = [i + 1/2) Ar with i = ■ ■ ■ — 1. The covariant spatial 
derivatives become covariant finite difference operators that are also associated with 
the links, e.g. 

exp[— ia,- Ar/2l </),• i 1 — <^,- , , 

Dr(f)^—^ — ^ — — i = 0---N-l. (39) 

where is short-hand notation for ai^j. 

It is now straightforward to discretize the action (|T3|) in a manner that still pos- 
sesses an exact local gauge invariance. But first, we need to state the restriction on 
4>o(t) corresponding to the boundary conditions ( |3lD and (pSj). Since oi is real, we 
can write these boundary conditions in a covariant fashion by requiring the real part 
of Dr(f) and the imaginary part of to vanish at r = 0. Using the discretized operator 



(39), we can then solve this boundary condition for 0o to obtain 

00 = Re[exp(-iao Ar/2)0i] , (40) 

where Oq is the value of ai i at z = and should not be confused with the time-like 
vector field. This now allows us to eliminate 0o from the list of dynamical variables. 
Finally, the discretized Lagrangian becomes 

^ I 2 V Ar J Ar2 J 

^"^ ' zao,j. , ,2 2 |exp(-zai Ar/2)0i+i - 0i|2 



+ 47r |l(5o-^ao,)x^l' + r,'l(5o- 
i=i ^ 

-I i\x.\' + 1)|0.P - M^X:<Pl) - ^(IX.P - 1)^ - Xr^ - l)^}Ar (41) 



2 [Im(exp(— i ao Ar/2)0i)]2 



-4^-1/2 , 

and the system may now be evolved using standard numerical techniques of ordinary 
differential equations. The Lagrangian ( ^Tj ) is actually of a Hamiltonian type with 
no dissipative terms, so it is convenient to use the leapfrog algorithm to perform the 
numerical integration. 

We do not have space to outline this well known computational procedure, so 
instead we simply state some of its more attractive features. First, the algorithm 
is second order accurate (i.e. the error from time discretization is of order (At)^ 
in the individual steps and of order (At)^ in an evolution of fixed length). Second, 
energy is exactly conserved in the linear regime, a desirable feature when pulling out 
the particle number. And finally, the algorithm possesses an exact discretized-time 
invariance, which is important since we are interested in obtaining the time reversed 
solutions starting from perturbations about the sphaleron. Of course these last two 
properties hold exactly only up to round-off errors, which can be made quite small 
by using double precision arithmetic. 



4 The Initial Configuration: Perturbation About 
the Sphaleron 



We are now ready to continue our investigation into the connection between the 
incident particle number of a classical solution and subsequent topology change. We 
could proceed by parameterizing linear incoming configurations of known particle 
number, but it would be extremely difficult to arrange the classical trajectory to 
traverse the sphaleron barrier. If we failed to see topology change for a given initial 
configuration, we could never be sure whether it was simply forbidden in principle by 
the choice of incident particle number, or simply because the initial trajectory was 
pointed the wrong direction in field space. 

To alleviate this difficulty, we have chosen to evolve initial configurations at or near 
the moment of topology change, and when the linear regime is reached the particle 
number will be extracted in the manner explained shortly. The physical process of 
interest is then the time reversed solution that starts in the linear regime with known 
particle number and subsequently proceeds over the sphaleron barrier. Of course we 
must explicitly check whether topology change in fact occurs, but we have found that 
it usually does. Fig. 2 illustrates the numerical evolution of the x fi^ld for a typical 
topology changing solution obtained in this manner. The modulus of x is represented 
by the height of the surface, while the phase is color coded (but unfortunately we 
can only reproduce the figure in gray scale). We have reverted to a gauge in which 
Xn = ~i and = 1, consistent with spatial compactification, and in which the 
incoming state has no winding and the outgoing state has unit winding number. The 
topology change is represented by the persistent strip of 2ti phase change near the 
origin after the transition. 

We turn now to parameterizing initial configurations. For classical solutions that 
dissipate in the past and future, topology change (and hence baryon number violation) 



is characterized by zeros of the Higgs field|ll]. For such topology changing solutions 



in the spherical Ansatz, the x field, which parameterizes the transverse gauge degrees 
of freedom, must also vanish at some point in its evolution. However, unless the 
transition proceeds directly through the sphaleron, the zeros of and x need not 
occur simultaneously, and for convenience we shall choose to parameterize the initial 
configuration at the time in which x vanishes for some nonzero r. Furthermore, we can 
exhaust the remaining gauge freedom by taking the initial x to be pure imaginary. We 
thus parameterize the initial conditions as an expansion in terms of some appropriate 
complete set with coefficients c^, consistent only with the boundary conditions and 
the requirement that x be pure imaginary with a zero at some r > 0. 

We choose to parameterize initial conditions in terms of perturbations about the 
sphaleron given by linear combinations of spherical Bessel functions consistent with 




Figure 2: Topology changing transition: behavior of the X field obtained the time reversal pro- 
cedure described in the text. The various shades of gray code the phase of the complex field. The 
field starts as an excitation about the trivial vacuum, passes over the sphaleron and then emerges 
as an excitation about the vacuum of unit winding. Note the persistent strip of 111 phase change 
near r = after the wave bounces off the origin. 

the small-r behavior (]26|)-(|3TD. We only need the first three functions 



since io{x) ~ 1, ji(x) ~ x and i2{x) ~ at small x. We also require the perturba- 
tions to vanish at r = L consistent with the large-r boundary conditions (|36| ) -(|38 D . We 
then parameterize perturbations about the sphaleron in terms oi jnmir) = jn{c(nmr) 
with n = 0, 1, 2, where anm with m = 1, 2, ■ ■ ■ are the zeros of jn{x). We are thus led 
to parameterize the initial conditions as 




(42) 



X 




(43) 



(44) 




(45) 



m=l 




(46) 



-^sph -^sph 

Xir,0) = 5Z C4m jlm(r) + Z 51 C5m J2m(r) (47) 
m=l m=l 

-^sph -^sph 

0(^,0) = '^(irn jomir) + i C7m jlmir) (48) 

m=l m=l 



sph 



aiir,0) = X! '^8m j2m(r) , (49) 



m=l 



where Xsph and 0sph are the sphaleron profiles, and where the sum is cut off at 
-^sph < N. To avoid exciting short wave length modes corresponding to lattice arti- 
facts, we shall take A'^^ph ~ (in our numerical work, A^sph = 50 for = 2239). 

We have used continuum notation, but (^)-(|i8|) is to be thought of as defining x 
and (p on the lattice sites and ai on the links rj+1/2. The time derivative of ai is to 
be determined by Gauss's law. 



5 Normal Modes and Particle Number 



We are now in a position to discuss the manner in which the asymptotic particle 
number is to be extracted. Recall that once the system has reached the linear regime 
it can be represented as a superposition of normal modes, and the particle number 
can be defined as the sum of the squares of the normal mode amplitudes. Since we 
have put the system on a lattice, we should properly calculate these amplitudes using 
the exact normal modes of the discrete system. However, since our lattice is very 
dense (A^ = 2239 with Ar = 0.04), it suffices to project onto the normal modes of the 
corresponding continuum system of finite extent L = NAr, the advantage being that 
we can solve for the continuum normal modes analytically. We have checked that this 
procedure agrees extremely well with projecting onto normal modes of the discrete 
system (obtained numerically), so for clarity we present only the continuum modes. 

It is convenient to work in terms of the gauge invariant variables of Ref . Jl3l . We 
write the fields x and (p in polar form. 



X 



-I [l + y] e'' 



(50) 



1 + 



h 



(51) 



where the variables y and h are gauge invariant. We can also define the gauge invariant 
angle 



(52) 



Finally, in 1+1 dimensions we can write 



r'f^, = -2t^,i) (53) 

where eoi = +1 and /i, v run over and 1, and where the variable ^l) is gauge invariant. 
Rather than working with the six gauge-variant degrees of freedom x, and we 
use the four gauge invariant variables p, a, and ^. 

We wish to find the equations of motion for small linearized fluctuations about 
the vacuum. In gauge invariant coordinates the vacuum takes the form h^^^ = y^ac = 
V'vac = ^vac = 0, and we thus need only work to linear order in the variables. From 
Ref. |T^ the normal mode equations are 




(54) 
(55) 
(56) 

. (57) 

Equation ([5^) corresponds to a pure Higgs excitation characterized by mass Mh = 
2\/X, while (^)- (|57|) correspond to three gauge modes of mass 

Note that there are four types of normal modes. The first two are easily obtained 
by solving the independent equations ( ^4|) and (H), while the last two can be found 
by solving the coupled equations (|56D and (|57D involving ip and ^. A solution in the 
linear regime can then be expanded as a combination of these four modes and the 
amplitudes akn, with k = 1,2,3,4 specifying the mode type, extracted. The Higgs 
and gauge particle numbers are defined by 

^ * mode 

^'higgs = X! 1*^1"! (^^) 

n=l 

N A 

* mode 



t'gauge 



{|a2„|^ + |a3„|^ + |a4„|^} , (59) 



n=l 



with total particle number given by 

y ^^higgs ~l~ '^gauge • (60) 

To avoid counting lattice artifacts we take the ultraviolet cutoff on the mode sums to 
be given by N^^dc ~ A^/5 to A^/10. 



^Upon restoring the factors of g and the Higgs vacuum expectation value w, these masses take 
the standard form Mh — \/2\v and Mw = 
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Figure 3: Decay of a small perturbation about the sphaleron: behavior of the particle number in the 
four modes as function of time for lattice parameters N = 2239, Ar = 0.04 and A^mode = 200 
with A = 0.1. The physical particle numbers are obtained by multiplying the asymptotic values in 
the graph by AlT / g'^ ~ 30, which gives A/higgs ~ 8 and A^gauge ~ 45, for a total physical particle 
number of A^phys ~ 53. 



Space does not permit a detailed exposition of this procedure, and one should 



consult Ref. |T0[ for full details. Here we must be content with Fig. 3, which displays 
the behavior of the particle number in the four modes of oscillation as a function of 
time. The initial state was a typical perturbation about the sphaleron as described 
in the previous section, and as it evolves it quickly linearizes and settles down into a 
definite asymptotic particle number. 



6 Stochastic Sampling of Initial Configurations 



Recall that our computational strategy consists in evolving a configuration near 
the top of the sphaleron barrier until it linearizes, at which point the particle num- 
ber can be extracted and the time reversed solution used to generate the topology 
changing process of interest. We can regard the energy e and the asymptotic particle 
number u as functions of the parameters c„ that specify the initial configuration, and 
by varying these coefficients we would like to explore the e-u plane and attempt to 
map the region of topology change. In particular, for a given energy e, we would like to 
find the minimum allowed particle number I'mini^) consistent with a change of topol- 
ogy. If this number can be made arbitrarily small, this would be a strong indication 
that baryon number violation would be observable in a two-particle collision. 

By randomly exploring the initial configuration space parameterized by the co- 
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Figure 4: Monte Carlo results with lattice parameters oi N — 2239, Ar = 0.04 (giving 
L=89.56), A^mode — 200 and A^sph — 50, and with a Higgs self-coupling of A = 0.1. The solid 
line marks the sphaleron energy egph = 47r(2.5426), below which no topology changing process 
can lie. The triangle represents the configuration from which wc seeded our Monte Carlo search. 
To obtain quantities in physical units, multiply the numbers along the axes by 47r/g^ 30. The 
energy axis extends from about 10 TeV to 15 TeV, while the particle number axis ranges from 
about 30 particles to 60. 



efficients c„, we would stand Httle chance of making lieadway. Instead, we sliall 
employ stochastic sampling techniques, which are ideal for tackling this type of multi- 
dimensional minimization. Our procedure will be to generate initial configurations 
weighted by W = exp(— F) with F = (3 e + /lu, and by adjusting the parameters /3 
and II we can explore selected regions in the e-u plane. In particular, by increasing n 
we can drive the system to lower and lower values of i' for a given e. In our numerical 
work we typically take j3 between 50 and 1000 while u ranges between 1000 to 20000. 

To generate the desired distributions wc have used a Metropolis Monte-Carlo 
algorithm. Starting from a definite configuration parameterized by c„, we perform 
an upgrade to c„ — > = c„ -|- Ac„ where Ac„ is Gaussian distributed with a mean 
of about 0.0008. We evolve the updated configuration until it linearizes and then 
calculate AF = (3 Ae + yU Az/. If the topology of the physically relevant time reversed 
solution does not change, then wc discard the updated configuration. Otherwise we 
accept it with conditional probability p = Min[l, exp(— AF)], which is equivalent to 
always accepting configurations that decrease F while accepting those that increase 
F with conditional probability exp(— AF). 

We are now in a position to present our numerical results. Fig. 4 represents 300 
CPU hours and involves 30000 solutions (of which only 3000 arc shown) obtained 
on the CM-5, a 64 node parallel supercomputer. We have chosen the lattice pa- 
rameters N — 2239, Ar = 0.04, with ultraviolet cutoffs determined by A^gph = 50 and 



-^mode = 200. The Higgs self-coupling was taken to be A = 0.1, which corresponds to 
a Higgs mass of Mh = 72 GeV. 

We have managed to produce a marked decrease of about 40% in the minimum 
particle number t^minle), which is approximated by the lower boundary in the Fig. 4. 
Nowhere, however, in the explored energy range does i/ drop below Att, or in physi- 
cal units the incident particle number N > 30 for energy E < 15 TeV (the outgoing 
particle number tends to be about 50 to 100). This is a far cry from two incoming par- 
ticles which would be necessary to argue that baryon number becomes unsuppressed 
in high energy collisions. 

The complex nature of the solution space can be illustrated by the break in popu- 
lation density between e/47r ~ 3 and e/Air ~ 3.4. In our first extended search we did 
not check whether topology change actually occurred, trusting that the time reversed 
solutions would continue over the sphaleron barrier. However, we later found an entire 
region between e/iir ~ 3 and e/47r ~ 3.4 in which the solutions never left the original 
topological sector. We excluded these points and restarted our search procedure near 
e/47r ~ 3. A small discontinuity in the lower boundary with slightly lower particle 
number was produced, but we have still managed to approximate ^'min(e) remarkably 
well. 

We can extract more information from the system by investigating the asymptotic 
spectral distribution |afc„p as a function of mode number n. Before we started the 
search, our seed configuration (represented by the triangle in Fig. 4) linearized into 
a distribution that was heavily peaked about a small mode number ripk ~ 50 (with 
An ~ 50), corresponding to a frequency of c<jpk ~ TTUp]^/ L ~ 0.1. After the search 
the solutions underwent a dramatic mode redistribution. The amplitudes |afc„ P of 
the linear regime peaked at higher mode number, ~ 75 — 100, with a much 
broader distribution (An ~ 200). Clearly our search procedure is very efficient in 
redistributing the mode population density. 

While V remains large throughout the energy range we have explored, it is in- 
teresting to note that t'min(e) maintains a slow but steady decrease with no sign of 
leveling off. To obtain an indication of the possible behavior of t'min(e) at higher ener- 
gies, we performed fits to our data using functional forms which incorporate expected 
analytical properties of the boundary of the domain of topology changing solutions. 
The fits gave a particle number iV = 2 at energies in the range of 100 TeV to 450 TeV. 
Of course we must explore higher energies before drawing define conclusions, but this 
is at least suggestive that particle number might at some point become small. 

While the energy range we have explored is of limited extent and more numerical 
work is clearly called for, it is still remarkable that we have extracted such a wealth 
of information from such an analytically intractable field theoretic system. Compu- 
tational techniques offer considerable promise in probing the nonlinear dynamics of 
the standard model, and we fully expect them to play a prominent role in our future 
understanding of high energy baryon number violation, both in extending the energy 
range of the classically allowed processes and obtaining information on the classically 



forbidden processes below the barrier. 
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